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Abstract 

We perform ab initio QED calculation of the (ls)^(2s)^2p3/2 - (ls)^(2s)^2pi/2 transition energy in the 
five-electron ion of argon. The calculation is carried out by perturbation theory starting with an effective 
screening potential approximation. Four different types of the screening potentials are considered. The 
rigorous QED calculations of the two lowest-order QED and electron-correlation effects are combined with 
approximate evaluations of the third- and higher-order electron-correlation contributions. The theoretical 
value for the wavelength obtained amounts to 441.261(70) (nm, air) and perfectly agrees with the experi- 
mental one, 441.2559(1) (nm, air). 

PACS numbers: 12.20.Ds, 31.30.Jv, 31.10.+Z 
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In the recent paper [ Ij, high precision measurements of the 2^3/2 — transition energy in Ar 
XV, Ar XrV, Ar XI, and Ar X have been reported. The best results have been achieved for the five- 
electron ion of argon, where the experimental precision is 2700 times better than the theoretical 
one and 200 times better than that in the best of the previous measurements. This unprecedent 
precision provides a unique possibility to test various branches of the theory describing many- 
electron systems. 

Since within the framework of the non-relativistic theory the 2^3/2 and 2pi/2 levels are degener- 
ate, the transition energy is exclusively determined by the relativistic and quantum electrodynamic 
(QED) effects. Hence, investigation of this transition allows us to test the many-electron QED ef- 
fects as well as calculations of the relativistic electron-correlation effects up to an extremely high 
level of accuracy. 

To date ab initio calculations of many-electron QED effects were considered for two- and three- 
electron ions 121131]. For systems with a larger number of electrons these effects were accounted for 
only within some one-electron or semi-empirical approximations |uJ,|4|,l5|]. Although the agreement 
of the results of Refs. Ill, |4i] with experiment is rather good, the precision of these calculations is 
much worse compared to the experimental one. The improvement of this precision is a challenging 
problem for the theory. The goal of the present letter is to improve the theoretical precision for 
the 2p3/2 — 2pi/2 transition energy in Ar^'^+. To achieve this goal we perform rigorous QED 
calculations to first two orders of perturbation theory and extremely large-scale configuration- 
interaction Dirac-Fock-Sturm (CI-DFS) calculations of the third- and higher-order contributions 
within the Breit approximation. 

To formulate the QED perturbation theory we use the two-time Green function method 
Instead of usual Furry picture, where only the nucleus is considered as a source of the external 
field, we have used an extended version of the Furry picture. It implies incorporation of some 
screening potential, which partly accounts for the interelectronic interaction, into the zeroth-order 
Hamiltonian. The perturbation theory is formulated in powers of the difference between the full 
QED interaction Hamiltonian and the screening potential. This accelerates the convergence of 
the perturbative series. In addition, the usage of the screening potential allows us to avoid the 
degeneracy of the (ls)^(2s)^2p3/2 and (ls)^(2pi/2)^2p3/2 states that occurs if the pure Coulomb 
potential is employed as the zeroth-order approximation. 

In the present letter we use four different types of the effective potential. The simplest choice 
is the core-Hartree (CH) potential. To obtain this potential we add the radial charge density distri- 
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bution of the (four) core electrons 

p, = 2Y,i9l + fn) (1) 

n=ls,2s 

with g and / being the upper and lower radial components of the one-electron Dirac wave function 
to the radial charge density distribution pnuc of the nucleus. The nuclear charge density is described 
by a Fermi distribution. The potential Vch generated by the total charge density p = pnuc + Pc is 
calculated self-consistently by solving the Dirac equation. 

The second choice is a local potential derived by inversion of the radial Dirac equation with 
the wave function obtained by solving the Dirac-Fock (DF) equation [Tl]. We will refer to this 
potential as local Dirac-Fock (LDF) potential. The construction of the potential Vldf is described 
in details in Ref. [|7|]. 

The other two potentials are based on the results of the density-functional theory (DFT). The 
first one is referred to as the Slater potential [SI]. This potential belongs to the wide family of Xa 
potentials. Introducing the total one-electron radial density pt via 

Pt(r) = 47rrV(r) , (2) 

POO 

p{r)cPr = I pt{r)dr = N , (3) 







where is the total number of the electrons, one can write the Xa potential in simple form: 



pt{r) a 



1/3 



V{r) = Kuc(r) + a ^ - ^rp,{r) . (4) 

Jo ^> r V327r2 j 

Here pt denotes the total one-electron density, i.e. includes both core-electron and valence-electron 

density, while the CH potential includes only the core-electron density. The value of the constant 

Xa for the case of the Slater potential is equal to 1 . To improve the asymptotic behavior of this 

potential at large distances, a self-interaction correction, known also as the Latter correction 

has been added to it. 

The fourth potential used in our calculations is known as Perdew-Zunger potential Vpz- It was 
constructed as described in Ref. [l^. 

The calculations of the transition energy can be conveniently divided in several steps. At first 
one has to solve the Dirac equation with the effective potential. Bound-state QED calculations 
require the representation of the quasi-complete set of the Dirac equation solutions. This was 
achieved by employing the dual-kinetic-balance (DKB) finite basis set method [lllll with basis 



functions constructed from B-splines 



Next we calculate the set of Feynman diagrams shown in Fig. 1 without any photon or electron 
loop, i.e. the part describing the interelectronic interaction. The dashed line ending with a triangle 
represents the interaction with the screening potential, taken with the opposite sign. The formulas 
for the calculations of the diagrams (a)-(d) in the framework of QED can be found in our previous 
works (see, e.g., Ref. [3]) devoted to the calculations of the two-photon exchange corrections to 
the energy levels of Li-like ions. We note that the diagrams containing only core electrons as initial 
(final) states can be omitted, because their contributions do not affect the transition energy. The 
formulas from Ref. [3] in our case should be completed by similar ones with the Is state being 
replaced by the 2s state. For the contributions of the diagrams (e)-(g) one can obtain: 

AE, = K., (5) 



AE, = 2 

c=ls,2s;ij,c 



where Kb = -(a|Kcr|&), hbcdi^o) = {ab\I{uj)\cd), I{uj) = e^af^a" D^^{uj), D is the photon 
propagator, P is the permutation operator, (—1)^ is the sign of the permutation, I'abcdi^) = 
{ab\^I{uj)\cd), and the indices c and v denote the wave functions of the core (Is or 2s) and 
the valence state, respectively. 

In the next step we should take into account the contribution of the diagrams depicted in Fig. 
[21 The evaluation of the one-electron QED corrections of order a, i.e. the self energy (SE) and the 
vacuum polarization (VP) (part (a) and (f) in Fig. ^ in an external Coulomb field is well known 



16h . High precision calculation of these diagrams for the case of an arbitrary external 



field, however, is a more involved problem y/Zi llSfl. We adopted the finite basis method for such 
calculations. 

The main QED contribution to the energy of the forbidden transition arises from the one- 
electron SE diagram. To reach the required precision, this contribution has to be evaluated with at 
least 0.1% of accuracy. With this purpose, we decompose the SE diagram into zero-, one-, two-, 
and many-potential terms, as indicated in Fig. [3l where the dashed line ended with a cross denotes 
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the full effective external potential. The zero- and one-potential terms have been calculated in 
momentum-space representation using the traditional renormalization scheme. The two-potential 
term has been evaluated in the coordinate space employing the analytical representation of the 
partial-wave decomposition of the free-electron Green function. To reach the required accuracy 
necessitates to sum over partial waves up to angular quantum number k = 50. Corresponding 
calculation within the B-spline approach would require an enormously large number of the basis 
functions, which would slow-down the computation considerably. To circumvent this problem we 
extracted the slowly-converging two-potential term and calculated it separately. A new numerical 
technique was elaborated for performing this calculation and will be described in details in our 
following work. The remaining many -potential term, containing three and more potentials, has 
been calculated within the DKB approach as the difference between the term, containing two- and 
more potentials and the two-potential one. The convergence of the partial- wave series for this term 
is very fast and, to reach the required accuracy, the summation can be restricted to \ k\ < 10. 

Since the contribution of the diagrams involving VP loops (parts (f-j) of Fig. |2l) was found to 
be very small, we calculated it within the Uehling approximation. The total VP contribution to the 
transition energy in the Uehling approximation is about 0.2 cm^^. The remaining Wichmann-KroU 
contribution is negligible. 

The other diagrams in Fig. |2](b-e) represent the self-energy screening contributions. The ir- 
reducible part of the diagrams (d) and (e) can be calculated as a wave function correction to the 
one-electron SE diagram. This part was calculated using the same algorithm as it was used for 
the calculations of the one-electron SE diagram. The diagrams (b) and (c) are known as vertex 
diagrams. They were calculated in the traditional way: The bound-electron propagator was de- 
composed into zero- and many -potential terms. The zero-potential term is ultraviolet divergent. It 
was renormalized and calculated in momentum space together with the zero-potential term of the 
reducible part of the diagrams (d) and (e). The remaining many -potential term contains infrared 
divergent terms. These divergences are canceled by the infrared divergences in the reducible part 
of the many-potential term of the diagrams (d) and (e) (see, e.g., Ref. UM)- AH the many -potential 
terms were calculated using the DKB approach. 

The contribution of the interelectronic-interaction diagrams of the third and higher orders has 
been evaluated within the Breit approximation using the CTDFS method. Going beyond the cal- 
culation performed in Ref. the basis set of the configuration state functions was significantly 
enlarged and the quadruple excitations were included. To separate out the contribution of the third- 
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and higher-order diagrams from the total energy, the following algorithm has been used. First of 
all, the Hamiltonian of the CI-DFS calculations was decomposed into parts 



H = Ho + XV, (8) 

Ho = -fffree + Kuc + Kicr , (9) 
V = V,,- \/,er , (10) 

where Hq is the unperturbed (zero-order) Hamiltonian, Kuc and Vee denote the operators of the 
electron-nucleus and electron-electron (Coulomb and Breit) interaction, respectively. V defines 
a perturbation and A is a freely varying parameter. This representation allows us to perform the 
expansion of the energy E in powers of A 

It is easy to see that the coefficients E^'^^ and E^'^^ correspond to the first- and second-order dia- 
grams depicted in Fig. [H calculated within the Breit approximation. Utilizing the same basis set 
these coefficients have been evaluated in two different ways: via numerical differentiation of E{X) 
with respect to A evaluated at A = and directly by means of perturbation theory. The term E^-^"* 
is then calculated as the difference between the total energy and the first three terms in Eq. (fTTI) 
evaluated at A = 1. The uncertainty of the E^-^^ contribution obtained in this way is estimated 
to be about 3 cm^. The QED contributions of the third and higher orders, which are beyond the 
Breit approximation, have not yet been evaluated. We estimate the uncertainty due to these effects 
as ±2 cm^^. 

The contribution of the one-electron two-loop QED-radiative corrections is very small and can 



be estimated using the analytical aZ-expansion reported in Ref. n2(M . Finally, one has to take into 
account the nuclear recoil effect. The calculation of this effect to all orders in aZ was performed 
in our recent paper [2]l|, where the isotope shift of the forbidden transition energies in B- and 
Be-like argon has been investigated. 

The results of our calculations are presented in Table 1. In the first line of the table the differ- 
ence between the one-electron Dirac energies of 2^3/2 and 2pi/2 states, calculated for the different 
screening potentials, is given. In the second row we give the contribution of first-order diagrams 
from Fig. [T] (diagrams (a) and (e)). These diagrams are calculated in the rigorous framework of 
QED, i.e. taking into account the energy dependence of the photon propagator. In the third line the 
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contribution of the diagrams of the second order from Fig. [T]calculated within the Breit approxima- 
tion is given. In the fourth line we give the QED correction to this contribution, i.e. the difference 
between these diagrams evaluated within the rigorous QED approach and within the Breit approx- 
imation. In the fifth line we give the contribution of the interelectronic-interaction diagrams of the 
third and higher orders derived from the CI-DFS calculations as described above. The contribution 
of the first- and second-order diagrams from Fig. [2] are presented in the sixth and seventh lines, 
respectively. The uncertainty due to uncalculated third- and higher-order QED effects is presented 
in eighth line. The following lines compile the contribution of the two-loop one-electron QED 
and the nuclear recoil correction, respectively. Finally, in the last line we present the total values 
of the energy of the forbidden 2^3/2 — 2pi/2 transition in B-like argon, calculated for the four 
different screening potentials. Averaging these values and accounting for the uncertainty due to 
the higher-order interelectronic-interaction and QED effects, we obtain i?tot =22656. 1(3. 6)cm^^. 
This value is four times more precise than that of Ref. fl: 22662(1 4)cm^^. The improvement 
mainly results from the calculation of the second-order QED effects (lines 4 and 7 in the table) as 
well as from the much more elaborated CI-DFS calculations performed in this work. Converting 
our result to the wavelength in air with the aid of Ref. [|22l] . one can obtain 441.261(70) (nm, air), 
which perfectly agrees with the experimental result 441.2559(1) (nm, air) from Ref. [1]. Further 
improvement of the theoretical value can be achieved by the reducing of the uncertainty of the 
CI-DFS calculations and computation of the QED diagrams of the third and higher orders. 

Summarizing, in this work we have calculated the energy of the forbidden 2^3/2 — 2pi/2 tran- 
sition in the five-electron ion of argon. The calculation incorporates the rigorous treatment of the 
second-order many-electron QED effects and the large-scale CI-DFS calculations of the third- and 
higher-order electron-correlation effects. This is the first ab initio treatment of the five-electron 
system in framework of QED. Significant improvement of the agreement between theoretical and 
experimental data has been achieved. 

This work was partly supported by RFBR grant No. 04-02-17574. A.N. A. acknowledges the 
support by INTAS YS grant No. 03-55-960 and by the "Dynasty" foundation. V.M.S., LIT., and 
G.R acknowledge the support by INTAS-GSI grant No. 06-1000012-8881. G.R also acknowl- 
edges support from DFG and GSI. 
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TABLE I: Various contributions to the energy of the forbidden transition in B-hke argon and the total result 
(in cm~^). 



VcH Vlbf Vpz Vsi 

Euirac 24343.0 25276.8 24860.7 26482.6 

E^^l -1002.9 -2064.2 -1539.3 -3451.8 

£'J2,Breit -1428.4 -1215.4 -1091.8 -993.3 

eSIqed 6.2 5.5 5.0 8.7 

^J^^Breit 694.0(3.0) 608.2(3.0) 376.0(3.0) 567.6(3.0) 

eJJ^d 47.9 49.9 48.9 52.5 

E^^l^ -2.2 -5.4 -3.8 -8.2 

E^^^l ±2.0 ±2.0 ±2.0 ±2.0 

£^two-loop -0.1(1) -0.1(1) -0.1(1) -0.1(1) 

^recoil -0.6 -0.6 -0.6 -0.6 

E^totai 22656.9 22654.8 22655.2 22657.4 
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FIG. 1: Interelectronic-interaction diagrams. 
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FIG. 2: Self-energy and vacuum-polarization diagrams. 
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